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We establish upper limits on branching fractions for B° decays to final states where the decay 
products are purely invisible (i.e., no observable final state particles) and for B° decays to vvj. 
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Within the Standard Model, these decays have branching fractions that are below current exper- 
imental sensitivity, but various models of physics beyond the Standard Model predict significant 
contributions from these channels. Using 88.5 million BB pairs collected at the T(4S) resonance by 
the BABAR experiment at the PEP-II e + e~ storage ring at the Stanford Linear Accelerator Center, 
we establish upper limits at the 90% confidence level of 22 x 10 -5 for the branching fraction of B° — > 
invisible and 4.7 x 10~ 5 for the branching fraction of B° — > vv^. 

PACS numbers: 13.20.He,12.15. Ji,12.60. Jv 



This paper describes a novel search for "disappearance 
decays" of B° mesons pj , where the B° decay contains no 
observable final state particles, or such "invisible" decay 
products plus a single photon. Invisible decay products 
are particles that are neither charged nor would gener- 
ate a signal in an electromagnetic calorimeter. These 
include neutrinos, as well as exotic, hypothetical parti- 
cles (such as neutralinos) . The rate for invisible B de- 
cays is negligibly small within the Standard Model (SM) 
of particle physics, but can be larger in several models 
of new physics. The SM decay B° — > vv, which would 
give such an invisible experimental signature, is strongly 
helicity-suppressed by a factor of order (m„/m B o) 2 0- 
When combined with the weak coupling constant Gp, 
the resulting branching fraction is necessarily well be- 
low the range of present experimental observability. The 
SM expectation for the B° — » vv^ branching fraction is 
predicted to be of order 10 -9 , with very little hadronic 
uncertainty [j| . An experimental observation of an invis- 
ible (+ 7) decay of a B° with current experimental sen- 
sitivity would thus be a clear sign of beyond-SM physics, 
as it could not be accommodated within SM theoreti- 
cal uncertainty. No quantitative experimental bounds on 
B° to invisible or 1/P7 have been previously established; 
however, a reinterpretation of data used for previous con- 
straints on b — > svv and other modes could potentially 
imply upper limits on the quark-level process of this de- 
cay 0. 

Several models of new physics can give significant 
branching fractions for invisible decays of the B° . A phe- 
nomenological model motivated by the observation of an 
anomalous number of dimuon events by the NuTeV ex- 
periment allows for an invisible B° decay to a Pxi final 
state, where x? is a neutralino, with a branching fraction 
111 the 10" 7 to 10" 6 range [falSl- Also, models with large 
extra dimensions, which would provide a possible solu- 
tion to the hierarchy problem, can also have the effect of 
producing significant, although small, rates for invisible 
B° decays 0,11 II 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II e+e~ collider. The data 
sample consists of an integrated luminosity of 81.5 fb _1 
accumulated at the T(45) resonance, containing (88.5 ± 
1.0) million BB pair events, and 9.6 fb" 1 accumulated 
at a center-of-mass (CM) energy about 30 MeV below 
BB threshold. The asymmetric energies of the PEP-II 
e + and e~ beams result in a Lorentz boost (3j m 0.55 of 



the BB pairs. 

A detailed de scrip tion of the BABAR detector is pre- 
sented in Ref. |lfj. Charged particle momenta are 
measured in a tracking system consisting of a 5-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer 
hexagonal-cell wire drift chamber (DCH). The SVT and 
DCH operate within a 1.5 T solenoidal field, and have 
a combined solid angle coverage in the CM frame of 
90.5%. Photons and long-lived neutral hadrons are de- 
tected and their energies are measured in a CsI(Tl) elec- 
tromagnetic calorimeter (EMC), which has a solid angle 
coverage in the CM frame of 90.9%. Muons are identified 
in the instrumented flux return (IFR), composed of re- 
sistive plate chambers and layers of iron that return the 
magnetic flux of the solenoid. A detector of internally 
reflected Cherenkov light (DIRC) is used for identifica- 
tion of charged kaons and pions. A GEANT4 based 
Monte Carlo simulation of the BABAR detector response 
was used to optimize the signal selection criteria and eval- 
uate the signal detection efficiency. 

The detection of invisible B decays uses the fact that 
B mesons are created in pairs, due to flavor conserva- 
tion in e + e~ interactions. If one B is reconstructed in 
an event, one can thus infer that another B has been 
produced. This technique has been exploited in several 
BABAR analyses 0, 0, [w| . We reconstruct events in 
which a B° decays to D^*^t + v, then look for consis- 
tency with an invisible decay of the other neutral B (no 
observable final state particles) in the rest of the event. 
In order to help reject non-BB background, R2, the ra- 
tio of the second and zeroth Fox- Wolfram moments, is 
required to be less than 0.5 flif . 

We reconstruct D^*'~ in the decay modes D~ — > 
K + it~it~ and D*~ — > D°tt~ where, in the latter case, D° 
is reconstructed in the decay modes K + ir~ , K + ir~ir°, or 
K + 7r~TT + 7r~ . To form D^*' candidates in these decay 
modes, K + candidates are combined with other tracks 
and 7T° candidates in the event. We identify K + candi- 
dates using Cherenkov-light information from the DIRC 
and energy-loss information (&E/<Ax) from the DCH and 
SVT. The 7T° candidates are composed of pairs of pho- 
tons in the EMC. Each photon must have a reconstructed 
energy above 30 MeV in the laboratory frame, and the 
sum of their energies must be greater than 200 MeV. The 
7r candidates must have an invariant mass between 115 
and 150 MeV/c 2 . A mass-constrained fit is imposed on 
7r candidates in order to improve the resolution on the 
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reconstructed invariant mass of the parent D meson. 

We require the D° and D~ candidates to have recon- 
structed invariant masses within 20 MeV/c 2 of their re- 
spective nominal masses except for D° decays with 
a 7T° daughter, which must be within 35 MeV/c 2 of the 
nominal D° mass. Mass-constrained fits are applied to 
D° and D~ candidates in order to improve the measure- 
ment of the momentum of each D. The difference in 
reconstructed mass between D*~ decay candidates and 
their D° daughters must be less than 150 MeV/c 2 . All 
candidates must have a total momentum between 
0.5 and 2.5GeV/c in the CM frame. 

Tracks selected as lepton candidates must pass either 
electron or muon selection criteria. We identify electron 
candidates using energy and cluster shape information 
from the EMC, and Cherenkov angle information from 
the DIRC. Muon candidates are identified using informa- 
tion from the IFR and EMC. Both electrons and muons 
must also have a momentum of at least 1 GeV/c, and a 
minimum of 20 DCH measurements. 

To select B° -> D^-g+u candidates, we require a 
£)W - candidate and a lepton candidate to be consistent 
with production at a common point in space. We then 
calculate the cosine of the angle between the £ + and 
the hypothesized B candidate, under the assumption 
that the only particle missing is a neutrino: 



cos 9 
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The energy E D (,)- e+ and mass m D (,)~ e + of the D^*'£ + 
combination are determined from reconstructed momen- 
tum information, and raj is the nominal B° mass. The 
B° momentum \pb\ and energy Eb are determined from 
beam parameters. When the assumption that a neutrino 
is the only missing particle is incorrect, cos &B,D(*)-e+ 
can fall outside the region [—1,1]. We thus require the 
£)(*}—£+ combination to satisfy —2.5 < cos6 B e>m-i+ < 
1.1. The asymmetric cut admits higher mass D* states 
where the additional decay products are lost, and al- 
lows for detector energy and momentum resolution. 
When more than one such D^~£ + candidate is recon- 
structed in an event, the one with the smallest value 
of | cos@ B £>(.*)-£+ 1 is taken. We reconstruct a total of 
126108 B° -> D^-£ + u candidate events in the data 
sample, with a purity of approximately 66%. 

We consider events with no charged tracks besides 
those of a B° — > D^*^£ + v candidate. Removing all de- 
cay products of the D^~£ + u candidate from consider- 
ation, we count the number of remaining EMC clusters 
consistent with a K® hypothesis, iV£ x tra , and with a pho- 

ton hypothesis, N° xtTa . Due to accelerator-induced back- 
ground and detector noise, the optimal requirements on 
iV£ x , tra and _/V° xtra are loose. For B° — > invisible candi- 
dates, we require that Np™ < 3 and N° xtIa < 3. For 
B° — » uDj candidates, we require only that there be one 
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FIG. 1: Distributions of U e xtra for (a) B° — > invisible and 
(b) B° — > uuj. The points with error bars correspond to 
data. The curves represent maximum likelihood fits to a sum 
of distributions modelling signal and background. 



remaining photon candidate with energy greater than 1.2 
GeV in the CM frame. 



The total energy in the EMC, in the CM frame, of pho- 
ton clusters that remain after the decay products of the 
D(*)~£+u candidate are removed, is denoted by -E e xtra- 
For B° — > vv^f, the energy of the highest-energy photon 
remaining in the event (the hypothesized signal photon) 
is also removed from -Eextra- In both B° —> invisible and 
B° — * uvj, this variable is strongly peaked near zero 
for signal, whereas for the background it is less strongly 
peaked, as seen in Fig. 1. The background can peak near 
zero due to events in which all charged and neutral par- 
ticles from the signal B° are either outside the fiducial 
volume of the detector, or are unreconstructed. For B° 
— * vv^j, the background shape increases at large E ext ra 
due to photons arising from misreconstructed ir° decays, 
and the best-fit amount of signal is slightly (but not sig- 
nificantly) negative. We construct probability density 
functions (PDFs) for the -E extra distribution for signal 
(.F s ig) and background (^bkgd) using detailed simulation 
of signal and background data. The background from ac- 
celerator and detector noise is modelled using randomly- 
triggered events in data. The two PDFs are combined 
into an extended maximum likelihood function £, defined 
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FIG. 2: Comparison of E cxt va between data (points with error 
bars) and Monte Carlo background simulation (histograms) 
for (a) B° -> invisible and (b) B° — > vv^j. The multiple 
categories of background in the detector (T(4S) — > B°B°; 
T(4S) B + B~; e + e~ — > cc; e + e~ — > uu, dd, or ss; and 
e + e~ — > t + t~ events) are each simulated by the Monte Carlo 
and plotted cumulatively. No background from e + e~ — > uu, 
dd, ss, or t + t~ is seen in the Monte Carlo sample. For 
both of the modes, signal would tend to peak strongly in the 
horizontally-shaded region. 



as a function of the free parameters iV s i g and iVbkgd 

e -(N sis +N bkgd ) 
£(N Big , iVbkgd) = ^ x 

N 

Y[(N siK F sie (Ei) + 7V b kgd^bkgd(^)), (2) 

where iV s i g and iVbkgd are the number of signal and back- 
ground events, respectively. The fixed parameters N and 
Ei are the total number of events in the data sample and 
the value of -E cx tra for the ith event, respectively. The 
negative log- likelihood (— ln£) is then minimized with 
respect to iV s i g and iVbkgd in the data sample. The result- 
ing fitted values of iV s i g and iVbkgd are 17 ± 9 and 19j^g° 
for B° -> invisible and — l.l±f"| and 28tf for B° -> vv^, 
where the errors are statistical. FigureGlshows the -Eextra 
distributions for B — > invisible and B a — > vv^. 

Using detailed Monte Carlo simulation of B° — > in- 
visible and vv^j events, we determine our signal effi- 
ciency to be (16.7 ± 1.0) x 10 -4 for B° — > invisible and 
(14.4 ± 1.0) x 10~ 4 for B° — > vv^j, where the errors are 
again statistical. For the B — > z/z/7 channel, we assume 
a photon momentum distribution predicted by the con- 
stituent quark model for B° — > z/z/7 decay, as given in 
Ref. Q. Of signal events that contain a reconstructed 
B° -> flW-fi/, approximately 46% (30%) of B a -> in- 
visible (_B° — * 1/P7) events pass the signal selection. 

We consider systematic uncertainties on the signal re- 
construction efficiency, and also the uncertainty on the 
ratio of background to signal determined in the fit. Sys- 
tematic uncertainties on the signal efficiency are domi- 
nated by the statistical size of the signal Monte Carlo 




FIG. 3: Comparison of -Ecxtra between data (points with error 
bars) and Monte Carlo background simulation (histograms) 
for the validation channels (a) B — * invisible and (b) B 
— > 1/1/7. No background from e + e~ — » t + t~ is seen in the 
Monte Carlo sample. 



sample (resulting in relative uncertainties of 6.5% and 
6.8% for B Q — ► invisible and B° — > wj> respectively) 
and by uncertainty on the efficiency for determining the 
particle type of charged tracks (5.4% for both channels). 
Systematic uncertainty on the number of signal events, 
due to uncertainty on the ratio of background to signal in 
the fit, is dominated by the parametrization of the back- 
ground and signal shapes (resulting in uncertainties on 
the number of signal events of 6.1 and 0.5 events for B° 
— > invisible and B° — > vvj, respectively) and by the en- 
ergy resolution for reconstructing neutral clusters in the 
EMC (3.2 and 3.4 events, respectively). Other systematic 
uncertainties include the efficiency for reconstructing the 
charged tracks in the B° — > D^~£ + v decay, the charged 
track momentum resolution, and the total number of BB 
events in the data sample. The total systematic uncer- 
tainties on the efficiency are 10.9% and 11.1%, and on 
the fitted number of signal events are 7.4 and 4.3 events, 
for B° —> invisible and B° — > 1/2/7, respectively. 

To determine 90% confidence level (C.L.) upper lim- 
its on the branching fractions of B° — > invisible and B° 
— > 2/2/7, we generate 8000 Monte Carlo experiments, each 
parametrized by the fitted numbers of signal and back- 
ground events, the efficiency, and the number of BB 
events in the data sample. Errors are incorporated into 
the simulated experiments via a convolution of the sys- 
tematic effects (treated as Gaussian distributions) and 
the statistical error (taken from the non-Gaussian likeli- 
hood function from the fit). 

The resulting upper limits on the branching fractions 
are 

B(B° -> invisible) < 22 x 10" 5 and 

B(B° -> 2/2/7) < 4.7 x 10~ 5 at 90% C.L. 

If the B° — > invisible branching fraction were zero, the 
probability of observing an equal or larger signal yield 
would be 6%. 
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We perform validation cross-checks on the results of 
this analysis. To check the measurement of the efficiency 
for reconstructing B° — » D^*'t + v decays (which was de- 
termined using Monte Carlo simulation) , we select a data 
sample in which a B° and a B° are both reconstructed 
as decays to D^tu in the same event. Using the ratio 
of such "double tag" data events to events where just a 
single D^*'iv is reconstructed, and the number of B°B° 
events in the full data sample, we determine the efficiency 
for B° — > D(*)~£ + is reconstruction in data. The result is 
consistent with that obtained from Monte Carlo simula- 
tion. 

We also search for the unphysical modes B^ — > invis- 
ible and B^ — > vDj (which would violate charge con- 
servation), to check that their resulting signal is con- 
sistent with zero. For these modes, we reconstruct B^ 
— > D°£vX°, where X° can be a photon, 7r°, or noth- 
ing. The D° is reconstructed in the same three de- 
cay modes as in B° — > D^*^£ + i>, and similar criteria 
are enforced for the reconstructed B as for the neu- 
tral B modes. All systematic errors are considered, and 
the "double tags" validation above is also performed 
for B^ reconstruction. The resulting fitted values of 
iV s ig are — 6^g°(stat.)±6(syst.) for B^ — > invisible and 
8i|(stat.)±4(syst.) for B ± — > vv'j, which are both con- 
sistent with zero. Figure shows the -£? e xtra distributions 
for the two validation modes. 

In summary, we obtain limits on branching fractions 
for B° decays to an invisible final state and for B° de- 
cays to 1/1/7. The upper limits at 90% confidence level are 
22 x 1CT 5 and 4.7 x 1CT 5 for the B° -> invisible and B° 
— > vv^( branching fractions, respectively. The latter limit 
assumes a photon momentum distribution predicted _by 
the constituent quark model for B° — > uvy decay |3j, 
whereas the B° — > invisible limit is not decay-model de- 
pendent. 
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